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Abstract 

Infrared growth of geometrical fluctuations in inflationary spacetimes is investigated. 
The problem of gauge- invariant characterization of growth of perturbations, which is of 
interest also in other spacetimes such as black holes, is addressed by studying evolution of 
the lengths of curves in the geometry. These may either connect freely falling "satellites," 
or wrap non-trivial cycles of geometries like the torus, and are also used in diffeomorphism- 
invariant constructions of two-point functions of field operators. For spacelike separations 
significantly exceeding the Hubble scale, no spacetime geodesic connects two events, but 
one may find geodesies constrained to lie within constant-time spatial slices. In inflation- 
ary geometries, metric perturbations produce significant and growing corrections to the 
lengths of such geodesies, as we show in both quantization on an inflating torus and in 
standard slow-roll inflation. These become large, signaling breakdown of a perturbative 
description of the geometry via such observables, and consistent with perturbative insta- 
bility of de Sitter space. In particular, we show that the geodesic distance on constant 
time slices during inflation becomes non-perturbative a few c-folds after a given scale 
has left the horizon, by distances ~ ~ RS, obstructing use of such geodesies in 

constructing IR-safe observables based on the spatial geometry. We briefly discuss other 
possible measures of such geometrical fluctuations. 
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1 Introduction 



A long-standing problem has been formulation of a more complete quantum treatment of in- 
flationary cosmology. Gravity has familiar short-distance problems {e.g. renormalizability, 
singularities), but also long-distance problems that become evident in inflationary cosmology, 
as well as in black holes. Those include apparent infrared (IR) growth of perturbations (see 
e.g. [l]), the general problem of specifying gauge-invariant observables (see e.g. [2]), and, in the 
former case, related problems of measure^ 

Understanding the IR growth of perturbations, particularly of the "gas" of gravitons pro- 
duced in inflation, has in part been confounded by the challenges of formulating observables 
sensitive to these perturbations. Interpretations of this growth has ranged from the claim that 
through backreaction they cause decay of the cosmologocal constant |4 -11 , to the belief that 



they have little meaningful effect. What is needed is a gauge-invariant way to describe physical 
effects of such fluctuations. 

Here it is important to note a distinction: if we imagine that there is a yet-unknown com- 
plete quantum-mechanical description of cosmology, then one of its basic features should be 
a set of quantum-mechanical observables. However, it is clear that, due to our limitations as 
Earth-bound observers in a particular era of cosmology, there is a much more restricted set of 
quantities that we can actually measure - and that an observational cosmologist would describe 
as "observable." Much of the discussion of this paper focusses on the former notion of observ- 



able. Where it is confusing, the former may be distinguished by calling them q-observables 12 
a given q-observable may or may not be observable to Earth's cosmologists. 

One goal of this paper is to formulate candidate q-observables sensitive to fluctuations in 
geometry. The latter in particular can have variance that becomes large at large times, and so a 
related goal is to better understand the meaning of this growth. The considerations of this paper 



are more formal than those of 12 , which discussed possible observable (to us) effects. These 
are related since an observable must ultimately be a q-observable. Moreover, there is common 
problem with both (related also to the measure problem) and that is to provide quantities that 



are "infrared safe," and thus well-defined even in the presence of large IR fluctuations. In 12 
infrared-safe observables were defined, relevant for the observations of a late-time observer of 
the CMB sky, by gauge transforming away metric fluctuations with wavelengths outside the 



horizon of the observer; a similar prescription was investigated in 13 



An alternate approach, in particular pursued in 16 20 , is that correlation functions written 



14 



in terms of the unperturbed geometric distance on the reheating surface would constitute 
^For a review see 3 . 

^In particular, in \L4\ it was similarly argued that one could apply such a gauge transformation, removing 
long wavelength modes. Their specific transformation leads to an exponential truncation of long wavelength 
modes given essentially by multiplying the IR part of the spectrum by a factor 1 — e"''*^ with p > 0. For arbitrary 
choice of p this does not correspond to the spectrum of the gauge invariant comoving curvature perturbation, in 



terms of which late-time observables are usually expressed by cosmologists 15 . However if p is identified with 
the size of the observed volume, then the proposal is essentially similar to the proposal of 12 , with the only 
difference that fl2] applied a step function instead of a smooth exponential window function. In the language 



of p corresponds to a renormalization scale (in 12 labelled l/g), set by the choice of observer 
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candidate IR-safe observables. This is not necessarily what a late-time observer measures, 
but is an interesting prescription to consider. Part of the problem with such a construction is 
that the curvature perturbation is not conserved when written in terms of the geodesic distance 
of the reheating surface. But even worse, as we show in section |6| when fluctuations of the 
geodesic distance are taken into account, this distance receives large perturbative corrections, 
and the description of the geometry in terms of geodesic distance becomes problematic. This 
happens only a few e-folds after a given scale has exited the horizon. 

In contrast, in terms of the comoving momentum the curvature perturbation is conserved on 
large scales, which makes the correlation functions of the curvature perturbation in comoving 
momentum at horizon crossing a convenient variable for discussing the observables for a late- 
time observer after the end of inflation. For this correlation function the effect of longer 
wavelength fluctuations, exiting the horizon even earlier, is to shift the background of the 
shorter wavelength mode when it exits the horizon. This can be described in terms of the semi- 



classical relations developed in 21 . To a late-time observer, who can measure the correlation 
function in many different patches of the sky, this becomes a physical effect. For a very late- 
time observer with access to a very large inflated volume set by of order N ~ l/H"^ e-folds 
(or a duration of order t ~ RdsSds, where Rds and Sds are the de Sitter radius and entropy 
respectively) of inflation, this effect becomes non-perturbative, while a present day observer 
with only access to the last 60 e-folds of inflation will see only a small effect on the CMB 
sky. A renormalization group equation relates the different observers with different survey 
volumes [12]. 

Despite the drawbacks from the viewpoint of present-day observation, q-observables based 
on spacetime geometry have some promise for giving gauge-invariant descriptions of inflation. 
For that reason, this paper will describe construction of q-observables of this kind, and investi- 
gate what they tell us about the physics of inflation. A specific focus will be on q-observables 
capable of probing growing variance, e.g. of tensor fluctuations, and interpreting its meaning 
for the geometry of inflation. In short, while we cannot support arguments for decay of the 
cosmological constant, rapid expansion produces buildup of long distance modes leading to 
large cumulative fluctuations after sufficiently long time. So, the growing variance does appear 
to indicate a quantum instabihty of inflationary spacetimes such as de Sitter (dS) space, to 
large excursions from the semiclassical geometry, and indicates a breakdown of our perturba- 
tive treatment. Indeed, this behavior is apparently another facet of the basic phenomenon 



responsible for self-reproducing inflation 21 



In outline, the next section describes in more detail the problem of observables in infla- 
tion, and the question of defining q-observables using geometrical measures such as lengths 
of geodesies. In particular, we note that the latter cannot be formulated based on spacetime 
geodesies spanning super-horizon distances, though one may base constructions on geometry 
of specific time slicings. We turn to study of these in subsequent sections. Section 3 sets up 
perturbative quantization of slow-roll inflation, in different gauge-fixings. Section 4 describes 
a basic construction, where one diagnoses properties of a fluctuating geometry by the lengths 
of curves in it. We note that growing metric fluctuations make important contributions, but 
also that care must be exercised in order to find gauge-invariant quantities. In particular, an 
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arbitrary curve is not gauge invariant, and as a consequence it is natural to consider instead 
curves that are geodesies of a fluctuating geometry. 

A simple context to examine such curves, and effects of metric fluctuations on them, is on a 
spatial torus T^, where we use a cycle of the torus to avoid the problem of specifying endpoints. 
Quantization on a toroidal version of dS is described. We find that tensor fluctuations do have 
a significant effect on torus cycles, giving a gauge-invariant statement of large corrections at 
large times, which we take to be a diagnostic of large fluctuations in the geometry of the 
spacetime. Section 6 then considers geodesies with endpoints, which may be thought of as 
quantum particles or classical "satellites," or as locations of field operators, in gauge-invariant 
constructions. We find similar large contributions of metric fluctuations, due to fluctuations of 
the geodesic in the geometry. These effects become particularly pronounced in slow-roll. This 
section also describes various subleading effects, quantum and otherwise, in such a construction. 
The large contributions of metric fluctuations to geodesic lengths signal breakdown of our 
perturbative description of the geometry via such q-observables. Section 7 considers the related 



problem of fluctuations in geometrical volume, making contact with 22 ,23 . Section 8 ends with 



brief description of other possible q-observable diagnostics of quantum inflationary geometry, 
and with conclusions. An appendix contains additonal details on gauge fixing. 



2 Q-observables and tensor fluctuations 
2.1 Challenges of gauge invar iance 

We seek gauge invariant observables in inflationary universes, particularly characterizing ge- 
ometry. Such gauge-invariant observables are surprisingly difficult to formulate in quantum 
gravity, where the symmetry of the effective field theory description is diffeomorphism invari- 
ance. This means that local observables familiar from field theory are not gauge invariant, and 
that observables must take a more non-local form. One would in particular like to formulate 
such observables that reduce to the local observables of field theory in an approximation; gen- 
eral aspects of this problem, and proposals for such "proto-local" observables, are described 
in|2). 

A focus of this paper will be on observables sensitive to fluctuations of geometry, and 
specifically to their apparent growth as evidenced e.g. by growth in variance of the tensor 
modes of the graviton. Here there are various puzzles. In particular, the growing variance 
is associated with tensor fluctuations that exit the horizon scale and are redshifted to long 
wavelength. However, at least locally a long-wavelength metric fluctuation is unobservable. 
Consider, for example, a small perturbation in a cosmological metric, with scale factor a{t), 

ds^ = -dt^ + a^{t){e^)ijdx'dx^ , (2.1) 

where is a spacetime-dependent metric perturbation (in notation of fT2]), which generally 
contains both scalar and traceless tensor parts. 

In the long- wavelength limit, Fjj is constant, and can be removed by a change of coordinates: 

x' ^ (e-^/2). .~i _ (2.2) 
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Figure 1: Metric perturbations can lead to significantly different geometries, as illustrated for 
the sphere. 



Indeed, this is reflected if one attempts to formulate observables in terms of local scalar quan- 
tities, such as the curvature; the curvature due to the perturbation vanishes in the limit of 



large redshift. This fact was also used in 12 , where it was argued that IR-safe quantities 



relevant to describing observations in a given observer's horizon can be formulated by using a 
transformation such as (2.2) to eliminate the effect of longer-wavelength fluctuations. 

However, there is a different viewpoint indicating that there is physics in long-wavelength 
perturbations of the gravitational field. In particular, the long-wavelength regime (compared 
to the Planck scale) is precisely the regime that gravitational wave detectors such as LIGO 
and Virgo are designed to probe. These detectors span kilometer sizes precisely because they 
are not measuring local quantities, but instead integrated quantities, specifically the relative 
displacement of two test masses separated by a finite distance. This physical picture suggests 
that more "global" observables are important in gravity. 

Indeed, consider spatial slices in an inflating spacetime; focus concretely on cases of closed 
slices, such as or T^. While it is true that the perturbations redshift as the slices expand, 
in order to compare two different slices one needs to perform a rescaling so that they have the 
same size. Once such a rescaling has been performed, growth of perturbations such as in (2.1) 
appear to lead to growing "lumpiness" of the space, as illustrated in fig. 1. 

There are various quantities one could try to use to characterize such growth of lumpiness. 
One possibility is to use an experimental apparatus to define an area, and measure rotation 
of vectors upon parallel transport around its perimeter; this kind of "Wilson loop" integral is 
sensitive to the curvature over the area's span. While the curvature may fall with the scale 
factor as 1/a^, the area of such a loop, if comoving, has opposite scaling, so the rotation of 
the vector remains fixed under scaling. One problem, however, is that specification of such a 
loop is not diffeomorphism-invariant, without additional physical structure. One could likewise 
describe other observable quantities (volume, etc.), though with similar problems of gauge 
invariance. 

An even more primitive kind of observable is the kind measured in gravitational- wave as- 
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tronomy, simply the distance between points defined by heavy point masses. These can be 
thought of as giving gauge-invariant "reference points" for measuring geometry. 



2.2 Q-observables via geodesic lengths 

We will pursue such an approach to defining observables sensitive to tensor fluctuations in 
cosmology. For example, one could consider a cosmology with a gas of massive test particles 
(or "satellites") present, and give one characterization of the gravitational fluctuations in terms 
of fluctuations of the distances between the test particles; these distances are given in terms of 
integrals of the metric over paths connecting the particles. 

A more fundamental characterization of this problem takes into account the quantum nature 
of the particles. Indeed, a related class of observables have been proposed, in order to reproduce 
the two-point function in field theory. Let a be a scalar field, and define the two-point function 
at a fixed geodesic separation S = s{x,x'): 

Such a construction has been investigated in (2| [T6[j20|[24] , and used in the context of the 



simplicial approach to quantum gravity 25 . This non-local construction can also be given for 
other local operators, such as the curvature scalar R{x). 

However, one immediately encounters an important limitation to characterizations of the 



geometry by lengths of geodesies, and of correlation functions via a construction like (2.3). 
This arises from the statement that arbitrary spacelike-separated points cannot be connected 
by a spacelike geodesic in inflationary geometries, and is most easily illustrated in de Sitter 
space. While we are familiar with the statement that in inflationary geometry, objects quickly 
become causally disconnected in the sense that there is no timelike or null geodesic connecting 
them, this property at first seems rather surprising in comparison to more familiar spacetimes. 
A quick argument is the following: suppose two points are connected by a spacelike geodesic. 
Then, by a dS symmetry transformation, the points and the geodesic may be taken to lie on 
the "equator" (minimal radius sphere) of the global slicing. But, no two points on this slice 
are separated by a distance greater than tt/H, so this serves as an upper limit on the shortest 
geodesic distance between two points. 



This property can be investigated in the flat slicing (2.1). Consider vanishing perturbation 
r, and, without loss of generality, a geodesic with tangent orthogonal to the = y and x^ = z 



directions. The F = metric (2.1) has Killing vector 

^^^ = (0,1,0,0) , (2.4) 
and so k^dX^/ds = c is conserved along a geodesic. Introducing the conformal time rj through 

dt = a{t)dri , (2.5) 



conservation immediately gives the equation 

drf = dx^ [l-'^] . (2.6) 



a2 



c2 
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For de Sitter space, 



and (2.6) can be rewritten 



a = e^* = -II {Hr^) , (2.7) 



dx = /^^^ , (2.^ 



with solution 



(x - xo)' = v^- vl ■ (2.9) 
This curve reaches maximum time a.t 1] = tjq < and x = Xq . Rewriting it as 

i/V(x - xo)' = 1 - 4 ' (2-10) 

we find that two such points in a time shce with rj < rjQ have a maximum "physical" distance 
measured in the slice given by 

aAx^ax = 2/i/ , (2.11) 
from the curve with qq ^ oo. The geodesic distance is easily computed, and found to be 

S = ^- ^tan-^ I 1 , (2.12) 

with maximum 5max = tt/H. Thus, for fixed comoving distance Ax, points rapidly exceed the 



bounds (2.11), (2.12), just outside the horizon, and there is no geodesic connecting them. This 
emphasizes the extreme causal and geometrical separation between points that is intrinsic to 
infiation. 

If there are no geodesies connecting generic points in an infiationary geometry, one must 
seek other ways to characterize the geometry. One approach which we will follow in this paper 
is to find geodesies subject to an additional constraint, e.g. that they lie in a given spatial 
slice, in a particular time slicing]^ This requires some additional structure, to play the role 
of a "clock." One candidate is an evolving scalar field, though without a potential an initial 
field velocity redshifts away. This does leave, however, the case of, e.g., slow-roll infiation. One 
can seek other candidates in the geometry itself, e.g. using the local expansion in the volume 
element as a clock. Or, one may consider additional physical structure. 



In a construction involving satellites or such as (2.3), one faces the question of disentangling 
the quantum dynamics of the particles from that of the gravitational field. For this reason, we 
will first focus on what appears to be an even simpler problem, temporarily dispensing with the 
particle "anchors" for the world-line ends, before returning to discuss these later in the paper. 
Namely, if we consider a classical de Sitter cosmology, in the fiat slicing, an equally good 
solution results from periodic identification under translations in to give the the torus T^. 
In such a geometry, the proper lengths of the cycles of the torus, measured via the integrated 



■^Another alternative is to consider curves with, e.g., fixed acceleration. We thank D. Marolf for discussions 
on such an approach. 
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distance, appear to be useful characteristics of the geometry. Passing to the quantum theory, 
we can consider gravitational fluctuations about our inflating torus, and see what their effect is 
on these length variables. Of course a proper treatment of this problem involves examination of 
possible gauge-dependencies of these length variables. However, we anticipate that this problem 
might be easier here, given that the curves are anchored to the homotopy of the torus, which 
is a topological invariant. 

Specifically, describe the torus via the identification ~ a;* + L. We will first consider as a 
candidate q-observable the length around the x-cycle of T^, given by the line integral 

S = <j) ds = <j) dx^Jgxx ■ (2-13) 

The zeroth-order in F result is S = a{t)L. However, fluctuations in Fjj will produce fluctuations 
in this length, and in particular we might anticipate significant corrections resulting from the 
infrared growth in these fluctuations. In order to explore this question, we need a careful 
treatment of their quantization. 

3 Action, perturbative expansion, quantization 
3.1 Action and perturbative expansion 

We begin with dynamical preliminaries. For concreteness we will focus on characterizing the 
geometry of models of single- field, slow- roll inflation, or of its "no-roU" de Sitter limit. The 
action is 



5" = 1 / d^Xy/^ 



^ (V0)2 - 2V{ct>) 



SttG 



(3.i: 



where TZ is the spacetime Ricci scalar. In the de Sitter limit, V is just the cosmological constant, 
V ^ K =const. 



A perturbative description of the coupled metric and matter fluctuations can be derived 27 
using the Arnowitt-Deser-Misner j28] (ADM) parameterization of the line element, 

ds^ = -N'^dt^ + hij{dx^ + N'dt){dx^ + N^dt) , (3.2) 

where N , are the lapse and the shift functions. 

We begin with a solution to the equations of motion, 

iV = 1 , AT* = , = a\t)5,j , = 0o(t) , (3.3) 

with instantaneous Hubble scale H = a/a, and satisfying the slow-roll conditions 

e = {y'f/2V'^ < 1 , V = y"/y < 1 • (3-4) 

An arbitrary perturbation about this may be parameterized as 

hij = a^(t)e^^(e^)ij , (p = (pg + (p , (3.5) 
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together with AN = N — 1 and iV*, where 'jij is a traceless matrix, 'ja = 0. (Here we adopt 
the convention of summing over lower indices with the flat-space metric 6ij.) 

One requires a gauge condition fixing the diffeomorphism symmetry; different choices will 
be discussed below. The lapse and shift act as Lagrange multipliers, whose equations of motion 
are the constraint part of Einstein's equations. After a gauge is fixed, one may solve these 
equations perturbatively to determine AA^ and iV* in terms of hij and 0. As a result, the 
dynamical degrees of freedom consist of the two polarization modes of the graviton, contained 
in 7jj, and one scalar, described by a gauge-dependent combination of ( and (j). In the de Sitter 
limit, this scalar decouples from the metric at quadratic order; alternately, with no field (p, the 
"time translation" symmetry implies that the scalar curvature perturbation ( can be gauged 
away. 

One treats the fluctuations in a perturbation expansion in 1/M2 = SttG; we will typically 
work in units where Mp = 1. In the slow-roll regime, ~ V/3Mp. 

3.2 Gauge fixing 

The gauge symmetries are diffeomorphisms, x'^ = x^' + e'^, under which the linear variation in 
the metric is 

= Vfe, + V(°)e^ . (3.6) 
This implies that the variations of 7^ -|- 2C5ij and take the form 

5(7., + = d,ej + d,e, + 2HeH,j + O {e\ £7, e<^, eN' , eAA^) 

6^ = e% + 0(e^e<^) , (3.7) 

where subleading terms are higher-order in the 1/Mp expansion. (Explicit formulas for these 
and higher order gauge-fixing are described in the appendix.) 

The traceless perturbation, 7jj, is independent of choice of time slice at this order. The 
spatial diffeomorphisms are typically fixed by choosing transverse gaugej^ 

da^, = . (3.8) 

Specifically, a spatial diffeomorphism gives the transformation (where we are careful to extract 
the traceless piece) 

d^l'ij = d^lij + die, + ^d,d,t, + O (e^ £7, eN' , eAA^) . (3.9) 
Thus, at linear order, we can satisfy 



by choosing 



da'i^ = (3.10) 

ei = — dk-fkj , (3.11) 



*Note that perturbations satisfying these conditions can be thought of as "gauge invariant" (see e.g. 29 ). 
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where is the Green function for = didi. One can use the expression (3.11) to fix the 

gauge (3.10) at higher order in 1/Mp as well, by working order-by-order and including the extra 
contributions in (3.9) on the right hand side of (3.11) (see appendix). 

Two commonly-used gauges for fixing the time-slicing are comoving gauge, 







(3.12) 



and traceless or spatially-flat gauge, 



C = o 



(3.13) 



As we see from (3.7), these choices can be made order-by-order in an expansion in 1/Mp. The 
choice (3.12) is very physical, in that constant-0 slices provide the time-slicing. The alternate 
choice (3.13) can also be characterized physically, as the choice where the evolution of the local 
volume element is given by the original unperturbed solution. 



Once the gauge is fixed by (3.10) and either (3.12) or (3.13), the lapse N and shift iV* are 
determined by the constraint equations. For example, in pure de Sitter, with the gauge (3.13), 
and at linear order, these are 



1 



1 



r(l) 







and 



3H 



which, when combined, imply 



{d'% + ^d,d,)N^'^ = . 



(1) 



(3.14) 

(3.15) 
(3.16) 



Thus, a linear-order solution is A^*^^-* = N-^'' = 0, and this procedure can be generalized to slow 
roll and higher-order contributions. 

Note also that the transverse-traceless choice (3.10) and (3.13) may not completely fix the 
gauge. In particular, we see that at linear level these conditions are unaffected by residual 
gauge transformations satisfying 



id's,, 



1 



didi)ei = 



(3.17) 



with a corresponding e'' chosen to preserve 7 = 0. The residual transformations are of course 
compatible with the constraints, as can be seen from the linear transformations 



(3.18) 



and comparing (3.16) and (3.17). The dS isometrics provide such residual transformations. 
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3.3 Quantization 



Once a gauge has been fixed, the action (3.1) can be expanded in the parameter H/Mp, 

S = So + S2 + Si. (3.19) 



Here, 5*0 is the action of the classical solution (3.3), and the first-order term vanishes by the 
equations of motion. The second order terms are given in, for example, 27 . For the tensors, 
we have in transverse gauge 



s. 



1 



2,t 



(3.20) 



The scalar action depends on the gauge. For example, in the gauge (|3.12|), one finds 



S- 



2,s 



2~ 



(3.21) 



Higher-order interaction terms described by Si can likewise be worked out, expanding in 
j + and in ip. In doing so, one solves the constraints for AA^ and A^*, order-by-order in terms 
of the other degrees of freedom; these then contribute to the interaction terms. 

Correlators are then computed perturbatively in these interactions, with the leading, gaus- 
sian, correlators determined by the actions (3.20 ), (3.21 ). In particular, these correlators exhibit 
a growing variance at long times; for example the two-point function for tensors, averaged over 
the proper length scale a/H to regulate the UV, behaves as [26] 



a{t) 



(3.22) 



where gives an IR cutoff corresponding to a "beginning" of inflation. 



4 Fluctuating line integrals 

It is well known that the variance of fluctuations of light fields in de Sitter can become large, 
but, as explained in the introduction, an important question regards the gauge-invariant or 
observable consequences of this statement. Naively, one might expect the inflated space to 
be very smooth. But, further reflection suggests that the accumulation of fluctuations at zero 
physical momentum, driven by the expansion, could have an important effect, connected to the 



growing variance (3.22). In particular, pile up of these fluctuations suggests that the spacetime 
could have a very inhomogeneous structure at the longest scales. 

We will study such inhomogeneities using the kinds of q-observable described in section 
[2} namely via integrated lengths of curves, which are basic characteristics of any geometry. 
Given the issues with unconstrained (four-dimensional) geodesies described in section |2| our 
focus will be on geodesies within the three-dimensional geometries corresponding to constant- 



time shces in either of the gauges (3.12), (3.13). Our basic picture is that space becomes 



10 



lumpy on large scales (see fig. 1), and the geodesies are sensitive to that lumpiness. We will 
investigate this question in two different contexts. The first is that where we imagine a pair 
of "satellites" (freely moving heavy masses), and we consider the distance between them. The 
second context, which is free of any such extra matter dynamics, is that where we imagine an 
inflating spacetime whose spatial slices have the topology of the torus, T^, and where we study 
curves with non-trivial homotopy on T^. 

In either case the basic expression for the distance is 

S{t) = / ^/h.jdX'dX^ , (4.1) 

integrated along a curve X^{t,a) in the spatial slice t =const. In the satellite case, Xj{t) are 
the locations of the satellites, and in the torus case we take Xi{t) = X2{t). The fluctuating 



spatial metric is given by (3.5). It is convenient to group the metric perturbations together by 
defining 

Tijix, t) = 7i,-(x, t) + 2C(x, t)6i, . (4.2) 

Also, in some cases it is useful to define a scale-dependent metric, as introduced in [l2] , 
which only incorporates fluctuations above a certain wavelength. To leading order in H, this 
can be defined by 

hij{q;x,t) = a'^(t)[exp{T{q,x)}]ij, (4.3) 
where Tij{q,t,x) is the scale-dependent metric fluctuation at scale q, 

Tij{q,t,x)= ——[2C{k;t,x)Sij + -fij{k;t,x)] (4.4) 



and ({k;t,x), jiji^k'jtjx) are mode functions with comoving momentum /cjj In what follows, 
we will assume q < a{t)H{t), such that this expression will be constant in time to a very good 
approximation. In [T2], a scale-dependent physical momentum was also introduced 



l^q,i{k, x) — [e ■'^ ', Hq (/c, x) — [c ''^'\ijkikj , (4-5) 

and it was then argued that leading IR-dependent higher-order effects are incorporated into the 
spectrum P by writing the tree-level two-point function Po^k) instead as a function of K,k{k, x): 

P{k, x)d^k = Po{Kkik, x))d^Kk . (4.6) 



The line integral (4.1) can be expanded in terms of the metric fluctuation Fjo-; if we use 



the scale-dependent metric (4.3), S also depends on the scale q determining which metric 
fluctuations are included. For example, suppose that we choose comoving coordinates on an 
initial slice so that the satellite separation or torus cycle are in the x direction. In that case. 



^At higher orders, the hmit in the integral (4.4) depends on F 12 
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we can compute the contribution of the fluctuations to the distance S along a line of constant 
y and z connecting the endpoints; this can be expanded as 



ds = a{t) / dx \J [e^] 



S = 

f I I 1 

= a{t) I dx{l + -r^jj, + -^xi^ix - -^xx^xx + ■■■) ■ (4.7) 

Here, for brevity we suppressed the arguments of Fjj. Next, we can compute the average effect 
of the fluctuations on this distance. In cases where our slices are chosen such that C = 0; then 
the linear contribution (F) vanishes]^ and we need to consider quantities quadratic in F to see 
an effect. (While this no longer holds in the case of slow roll in a slicing with C 7^ 0? later in 
the paper we will demonstrate slow-roll suppression of such contributions.) One possibility is 
to look at the variance of 5, via 

{SS) . (4.8) 

The contribution to this from the linear term in F can be computed, but is found to be small 
(see section 6.2.4). Intuitively, the linear contributions to S average out when integrated around 
the x-cycle. 

The fluctuations also enter at quadratic order, suggesting a growing contribution to {S). 
This can be written in terms of the variance of the fluctuations as 

{S) = a(t)L (l + i(C^(a:)) + \{^^,{xHx{x)) - + • • • ) • (4.9) 



Eq. (3.22) and the analogous equation for C, show that these contributions grow with time. 
However, while the basic message of this equation - growing contribution due to the variance - 
is correct, this equation is not particularly physical as it stands. The problem is that the curve 
y =const., z =const. is rather arbitrary from the viewpoint of the perturbed metric. This has 



the effect of making the result (4.7) gauge-dependent, hence unphysical. 



Specifically, we could have described this computation by still working at constant y and z. 



but in a different gauge from (3.10) and get a different result. The general gauge transformation 



corresponds to the action of the diffeomorphism x^ — )■ x^' = x^ — e^{x) on the metric. However, 



from the diffeomorphism symmetry of the basic expression (4.7), we can alternately describe 
the situation by working with fixed metric, and consider a deformation of the curve; this 
corresponds to working in the x' coordinates, and has the effect — )■ x'^ + on the curvej^ 

To make a gauge- invariant statement, we would like a quantity independent of such a 
variation, which we rewrite in general as 

X'^(a) = X^{a) + 5X^((t) , (4.10) 
^Validity of such gauge fixing to second order is checked in the appendix. 

^What we have described is the standard switch between active and passive viewpoints, namely from a 
transformation on the fields to a transformation on the underlying coordinates. Note also that in the case of 
the path connecting satellites, the diffeomorphism also acts on the satellite trajectories. 



12 



where Xq (a) is the original path. In our particular example 

Xo^ = (0,x,0,0) . (4.11) 

We will first focus on the purely spatial gauge transformations (5X*, viewing e° as fixed by 
choice of physical "clock," determined by either of the conditions (3.12) or (3.13); we will 
examine these in particular cases. While introduces arbitrariness into 5, one might expect 
that there is a unique choice, corresponding to a geodesic in the new metric. Then, one could 
calculate the average variation of the length of this geodesic, as a gauge invariant. 

In order to explore this, we consider the change of S under a combined variation and 



hij — ho_ij -\- dh 



0,ij 



(4.12) 



where /lo.ii is given in (3.3) (we work at fixed time). Expanding to second order, this change 
takes the form 



S[Xq + 6X, h + Sh] 



+ 



5S^^ 5S 

Th'^^Jx''' 



1 



6X^ + 



6^S 



5X5h + ... 



(4.13) 



2 6h^ 2 5X'^ " ' 5X5h' 
where integrals are implicit and all variational derivatives are evaluated at Xq, h^. The linear 
term in 5X vanishes because the original curve was a geodesic. The linear and second-order 
terms in 5h combined are given in terms of and C, by (4.7). Our problem, therefore, is to 
understand the remaining 5X-dependent terms. 

Following the preceding discussion, we can try to fix 5X by the condition that the perturbed 
curve is a geodesic in the perturbed metric. Working to linear order in both perturbations, this 
is the condition 







5S 

Jx 



Xq+SX, 
hQ-'rSh 



5X2 



6X' + 



6X6h 



5XSh 



with second derivatives evaluated at Xq, /iq- 

.^±Plx 

SX dsl 



This can be rewritten 
5 D'^Xq 



-6h 



6h ds"^ 



(4.14) 



(4.15) 



found by varying the usual geodesic equation. For a solution of (4.14), (4.15), the last two 
terms of (4.13) may then be written in terms of 



5X 



:5X' 



ds 



d6X' d6X^ 
ds ds 



rfX, 



ds 



dXl 
ds 



For our x-path (|4.11|), the geodesic equation becomes 

1 



■pa 



(t) ( -dJKa + ^d„6h. 



a 



y,z 



(4.16) 



(4.17) 



where in the latter equation we find the linearized Christoffel symbols. A deformation 6X^ 
corresponds to a reparameterization of the path and is unconstrained. We will study the 
problem of solving these equations both for and for the satellites, and of determining the 
resulting change in path length in (4.13). 
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5 Cycles on the torus 



5.1 Quantization on the torus 

For simplicity we begin with the case of the spatial manifold T^, and with no matter present. 
We take the background metric to be (3.3), with the scale factor of dS, (2.7), now with the 
identifications ~ + L. The preceding discussion of fluctuations and gauge fixing thus 
carries over with minor modification|3 

To investigate the possibility of residual gauge transformations, consider the positive definite 
integral 



< d^x 



{die j) {die j) + -{die,i){dj 



e,\ 



(5.1) 



where the last equality arises from integration by parts, with no surface term for single-valued 
e*. For a residual gauge transform, the latter integral vanishes by (3.17), and then by positive 
definiteness we find that die^ = 0. Thus, the only residual gauge transformation is a trivial 
shift by a constant. 

In the transverse-traceless gauge (3.10), (3.13), the quadratic action is simply that for 
tensors, (3.20). The corresponding general linearized solution is given by the mode expansion, 

1 



E Ez^ K<,i^hk{t)+b^\et;{-khm 



(5.2) 



Here the first term is a space-independent zero mode contribution with mode index /. The 
momentum sum is over quantized values k = {27i/L){nx,ny,nz), with integer nf, specifically, 
quantization is related to that in infinite space by the replacements 



(27r)3 ^ Z3 



{27i)^6^ (k - k') ^ L^6] 



kk' 



(5.3) 



The ladder operators with s denoting the helicity, therefore satisfy the commutation relations 



"k' "k' 



(5.4) 



The polarization tensors are chosen to satisfy the transversality and tracelessness conditions 
e|j(k) = kie^{\^) = 0, together with a completeness relation obtained by tracing over spatial 
indices, efj(k)e*j'(k) = 25ss'- The mode functions are the same as those for a massless scalar, 
up to normalization; they are most easily written in terms of the conformal time 



Ha{t) 



as 



lk{v) 



H 



[1 + ikri)e' 



-ikrj 



(5.5) 
(5.6) 



^For previous discussion of quantization on T'^, see jsj. 
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The quantization of the zero modes is handled separately. In the case of finite volume, these 
modes have finite norm. The equation of motion for them has spatially homogenous solutions 
of the form 

70/ = Qi + Pi [^^ (5.7) 



with conjugate variable 



{5.t 



In particular, the canonical commutation relations for 7, together with completeness, imply 
commutation relations of the form 

i 

V 

where V is the volume of the compact space, and we drop the mode indices. The "dS-invariant" 
zero- momentum vacuum P\P = 0) = is non-normalizable just like the ground state for an 



[Q,P] 



(5.9) 



ordinary free particle 30 ; the corresponding wavefunction is 70 independent, so it gives a 



completely indefinite 70. 

To specify an initial size/shape for the torus, we instead choose the vacuum wave function 
for the zero modes 70 to be gaussian wave packets peaked around some definite values (here 
taken to be zero) in position space 



27r(AQ) 



1/4 

2 , exp 



(5.10) 



The initial velocity due to the finite width of the wave function quickly redshift away, and 
since the energy density of the zero mode redshifts away like a^7Q oc it can effectively be 
ignored. 

This type of initial state belongs to a class of 0(4) invariant Hadamard Fock vacuum 



states 31 , which can be constructed by defining 

p = (^Aao + A*al^ 



(5.11) 



where the commuation relation, 



1, implies 



A*B - B*A 



The Hadamard Fock vacua are then defined by 



ao|0) 



HF 



V 



. 



(5.12) 



(5.13) 



With the choice A = -i/{y/2VAQ) and B = AQ/V2, the vacuum wave function in position 



space reduces exactly to the one in eq. (5.10); see 31 
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5.2 Growth of variance 



The two-point function gives an important measure of the gravitational fluctuations; specifically, 
as in (3.22), consider the double trace. 



(5.14) 



Here, the IR divergence of the analogous expression for non-compact space has been regulated 
by the finite size of the torus, which plays the role of a comoving IR cutoff, like was used 
in [l2||2l||32]. The zero-mode contribution asymptotes to a constant oc (AQ)^ at long time. 
The coincident limit x = x', which gives the variance, is UV divergent, but may be regulated 
by choosing a minimum physical separation a(t)|x — x'| ~ ^/H, effectively providing a UV 
cutoff at ~ 27ia{t)H. This introduces a logarithmic time-dependence into the variance, as in 
the non-compact case; specifically. 



2(- 

V27r 



2na(t)H 



1 + A;V) 



2tt/L 



2^1og[a(t)ifL] 



(5.15) 



which grows linearly with time. We will also need the variance of specific tensor components. 



which can be derived using the sum over polarizations 33 



^ij,kl\^) 



Sijhi (5.16) 
+ SijQkqi + ^kiqiQj - SikQjqi - Suqjqk - SjkQiqi - Sjiqiqk + mjqkqi 



where q is the unit vector in the direction q. For an arbitrary unit vector n, we find 21 

4 8 
niUji-fiaij) ^ -(7^ (a;)) , runjUknii-fij-fki) ^ —(7^(2;)) . 



(5.17) 



We can use (5.17) to evaluate the tensor variance in (4.9); for generality we likewise include 
the scalar contribution, although in the present section it is gauged to zero. The result is that 
on the unperturbed path the line integral yields 



{S)=ait)L{l + l{eix)) + ^{^'ix)) + ... 



(5.18) 



In de Sitter space or its version, the variance, (7^(x)), keeps growing linearly with time 
like HH at leading order, and the effect becomes order one at a time scale of order t ~ RdsSds- 



(This is in parallel with the apparently related case of black hole evolution [21,34,35].) But 
even at the end of the non-eternal regime of chaotic inflation, we have (7^(x)) ~ 1 at the time 
of reheating, if initial conditions are set at the end of the self-reproduction regime. 
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Figure 2: Pictured is a torus, made "lumpy" by metric perturbations, along with a representa- 
tive cycle along which distance S can be calculated. 



5.3 The line integral on the torus 



As argued in section 4, the line integral along the unperturbed path does not give gauge- 
invariant information, and it is no longer a geodesic in the perturbed metric. We can try to 
remedy these by trying to solve the geodesic equation in the constant-time hypersurface to 
second order, in order to find the curve with the shortest distance around the torus. 



The geodesic equation (4.17) can be integrated, giving 



drSX' 



a — 



^ XX 1 



(5.19) 



with integration constants a". However, for a closed geodesic, (5.19) should be periodic, which 
is the condition 



From the expression (4.17) we find that this requires 



da 



dx Shr, 







(5.20) 



(5.21) 



Let us first consider the case with vanishing fiuctuations 6hxx = 0. Then a perturbed 
geodesic can be found by integrating and enforcing periodicity. It is useful to introduce a 
notation for an average over the curve, 

1 



[/] 



L 



dxf{x) 



(5.22) 
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with f the integral from to L, in terms of which we find 



= 5X^ = 



(5.23) 



up to an overall constant shift. We can then evaluate the total effect on {S); since the geodesic 
equation is satisfied and can be written as (4.14), we find that the last two terms of (4.13) have 
average 



I 5^S 



\5XSh 



6X6h) 



6^S 
2 \6X^ 



6X' 



a{t) 



dx {dJXadJXo,) (5.24) 



[c.f (4.16)). Using (5.23) then gives 



2 \5X2 



2a3(t) 



(5.25) 



The last term is subdominant, due to cancellation of the fluctuations in the path average, as 
in ( |48| ). 

Next, let us find the additional terms present when 5hxx 7^ 0. First, there are the quadratic 
terms in 5h already present in (4.9). Then, there is also a cross term linear in 5X and 

1 



2a{t) 



dxda.5hxx5X^ 



(5.26) 



If we write 5X = 5Xi + 6X2, with SXi given in (5.23), then the term involving 5Xi averages 
to zero since {'^xx'^xa) = for a 7^ x. There are also quadratic terms in 5X. If 5Xi and 8X2 
are uncorrelated, {dx^X^dx^X!^) = 0, so we have a term as in (5.24) involving only 6X2. 

Thus, combining these terms, the first, leading, term of (5.25) and the other terms in (4.13) 
gives 



(S) = a{t)L 



1 - 



dx 



1 



2a{t) 



^ {ixilix) - ^ {ixxlxx)^ + O ((7^(3^))^) 

ait) 



dxdaShxxSXn ) + 



dx {dx6X°dx6X^) 



(5.27) 



Consider the case 6X2 = 0. Using eq. (|5.17|), we find 

{S) = a{t)L 
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^\x)) + O U^\x)) 



(5.2^ 



in the spatially-flat time slicing (3.13), with (7 ) as in (5.15). Comparing with eq. (5.18), 
we observe that, while the circumference around the original curve becomes larger due to the 
growing variance, we can always find a nearby curve that is shorter than the original curve, by 
an amount that grows with the variance of the metric. 



18 



For general 6hxx there will not be a neighboring geodesic, but the length of the curve can 
be reduced further by choice of 6X2. For example, with constant 6X2, S changes by 

^ 6X^ (f dxdjh^^ . (5.29) 



2a{t) 

So, when da § dxShxx 7^ 0, {S) is further reduced by 8X2 with a particular correlated sign. 

An intuitive picture for our statements is the following. First, variations Sh^y, Sh^z cor- 
respond to a contraction of space in some directions, and expansion in other directions. The 
basic idea is that by choice of the path, we can take advantage of the contraction and shorten 
the total path length. On the other hand, consider a variation 6hxx{y, z). This has the effect of 
changing the circumference of the x-cycle in a ?/ and z dependent fashion, as sketched in fig. 2. 
We can decrease the length of the curve by sliding it towards smaller circumference, although 
reaching a minimum of the circumference, corresponding to a geodesic, may require a large 
deformation. 



But, since the length can always be lowered below (5.28), we have an invariant statement: 



for some curves in the relevant homotopy class, the average {S) is reduced by at least the 



amount in (5.28), and this reduction grows with the variance of 7. 



6 Satellite distances and correlator modifications 

An alternative way to make diffeomorphism invariant statements is to consider a length between 
two physical comoving particles in an inflationary spacetime; in the classical limit we might 
also describe them as "satellites." We assume that at some initial time they are separated by 
comoving distance L, in the x direction. 

If we compute the distance between the satellites as a function of time, a number of effects 
besides the inflationary expansion can contribute. First, the satellites can have an initial 
velocity. Classically, this may be set to zero, but in the quantum case it has nonzero spread 
due to position/momentum uncertainty. Secondly, with zero such velocity the satelhtes in 
the zeroth-order metric stay at fixed comoving location in the (x, z) coordinates, but metric 
fluctuations drive motion of the satellites. Thirdly, at sufficiently long times new "Boltzmann" 
satellites that are indistinguishable from the original satellites can be created, confounding the 



definition of the observables 36 . Finally, we have contributions of the metric fluctuations to 
the path length between the satellites, both due to local changes in the length, and due to 
the need to perturb to a new path satisfying the geodesic equation in the new metric on a 
constant-time slice, as described in the preceding sections. 

As we will argue in subsection |6.2 the inflationary expansion rapidly redshifts away an 



initial velocity, as well as gradients responsible for satellite motion, and so these effects, as 
well as those of Boltzmann satellites, are small on timescales where metric fluctuations make 
significant contributions to the path length. For that reason, we will first consider the latter 
set of effects. 
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6.1 Perturbations in path length 



We begin by investigating the change in length of the path between two satelhtes with fixed 
comoving coordinates. Thus, we want the contributions of metric and path fluctuations to the 
integral (4.1), with endpoints taken to lie at Xi{t) = (0, 0, 0) and X2{t) = {L, 0, 0). In the case 



of slow-roll with the comoving gauge (3.12), there is a contribution due to {() ^ 0; however. 



we will see in the next subsection that this is small relative to the other effects we consider. 
To second-order in fluctuations of the metric and path, the change in length is again given 



by (4.13). In the satellite case, the geodesic equation (4.14) can be solved for the geodesic 



connecting the endpoints. This allows us to eliminate the cross term proportional to 5h5X in 



(4.13). The result for the average {S) is 



(5) = 5o + a{t) 



dx 



1 



^ xiJ^ ix) Q xx^ xx") 1 



a{t) 



dx (dxSXadxSXa) + 



(6.1) 



4 8 / 2 

In order to evaluate the last term in this equation, we must solve the geodesic equation, which 



to linear order in perturbations is given by (4.17). The solution is of the form 



= - / dx' / dx'T'^^ + + 
10 Jo 



(6.2) 



where again a = {y,z). The boundary condition at x = gives 6" = 0, and the boundary 
condition ai x = L flxes a". If we introduce the notation 



r = - 



dxT' 



and use the path-average notation (5.22), the solution becomes 

6X" = I dx'ir - in) ■ 



(6.3) 



(6.4) 



Thus, the second-order variation in the path length, (6.1), is fully determined in terms of the 
two-point functions of the metric fluctuations. 

The fact that ( and 7 are independent random variables at quadratic level is sufficient for 
us to be able to minimize the path in ( and 7 independently. Let us therefore consider the 
physical time-slicing given by = (po, (3.12), and first focus on the effect of 7 fluctuations are 



subdominant in the slow- roll expansion. In this case we have leading contribution (F^jFja.) 
(r..F,,) =4(C2), thus 



(5) =So + lait) r dx {C) Tdx {dJXMX^) . 

^ Jo Jo 



(6.5) 



Let us flrst consider the situation where we neglect the perturbation of the path, as was 



done in |l7|[l8|[20j (see e.g. sect. IV.B. of [17) or eqs. (66), (67) of [18]) and [19]. Then we 
simply have 



{S) = a{t)L ( 1 + i (C^> + 



(6.6) 
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As described in section [2| we might consider the two-point correlation function, as a function 
of path length. In terms of coordinate separation, this correlation function is 



(C(o)C(x)) 



(27r)3 4^3 



n»-l 



(6.7) 



If we take x = (L,0,0), and write the correlation function in terms of (S) instead, then we 
obtain 



ris-l 



where we defined k = ke In 17 19 this object was considered as a candidate IR-safe 

observable, although this construction is not necessarily relevant for a late-time observer of the 
CMB sky for which an IR-safe observable was instead constructed in 12 . One problem with 



the construction above from the point of view of a late-time observer is that the curvature 
perturbation is not conserved on super-horizon scales when written in terms of the geodesic 
distance, since the geodesic distance receives large corrections on large scales. This implies 
that in terms of the of the average of the actual distance (S) along the original path on the 
reheating surface, neglecting the perturbation of the path, the spectrum become^ 



Pdk) = P^'\k) 



l + -{ns-l){C) + ... 



(6.9) 



where only modes that exit the horizon after the mode k, but before reheating, contributes to 
the variance. Thus, for a cosmologically relevant scale exiting the horizon only 60 e-folds before 
the end of inflation this effect is order H^N, which is small. But for very large scales, which 
exited the horizon ~ e-folds (or At ~ RS) before the end of inflation, the effect is 

large 
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Moreover, we have argued that the original path does not have any invariant meaning in 
the perturbed metric; the natural choice is to perturb to the new path which is a geodesic in 
the perturbed metric. We will again find this has a significant effect. 

Specifically, again focus just on the dominant scalar fluctuations, with 



■pa 

XX 



Using the notation (6.3), we now have 



fa{x) = I dx'daC{x') 



(6.10) 



(6.11) 



^The same formula was given by [37[. 

-•^^Note that this is different than the effect studied in |2f . There new semi-classical relations was used to find 
the effect of long wavelength modes on the correlation functions written in comoving coordinates (not in terms 
of geodesic distance on the reheating surface). In comoving coordinates the curvature perturbation is conserved 
on super horizon scales, but long wavelength modes, exiting the horizon before the observed mode, shift the 
background of the observed mode when it exits the horizon as in (4.6). To a late-time observer, comparing 



correlation functions at different scales, this becomes a physically observable effect 12 
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Thus, the contribution in (6.5) from the path perturbation becomes 



a{t) 



dx {dJX^dJX^) 



a{t) 



(6.12) 



We then insert the mode expansion for C,, and obtain 



d^k kl 



(1 



Next, we integrate over x, and define a variable r = kx/k, so that A;, 
give 



(CkC-k) 



k^(l-r^ 



(6.13) 
These steps 



dx (ifaix) - [fa]f) 



k'^dk 



(2 



TT 



(CkC-k) 



dr 



1 - {kLrf + 2cos{kLr) - 2 
k'^r'^L 



(6.14) 



At large /c, the integral over r in (6.14) behaves as ~ nkL'^/G. One can readily see that 
the dominant contribution to this behavior comes from the regime of small r, that is k^ k, 
\ka\ ~ k. For spectral index ng > 0, this produces a UV divergence in (6.14). For example, in 
the scale-invariant case Hs = 1, with (CkC-k) = ^[1 + {k/aHY]/k^, and introducing a comoving 
UV cutoff Auv, the contribution to the path length becomes 



a(t) 



dx{iUx)-[U]f) = a{t)L 



A 

48^ 



Ar/y + 



3{aHy 



L + 



(6.15) 



where the dots represent terms that grow less rapidly with the cutoff. 

Note that working with the scale-dependent metric (4.4) with scale q corresponds to taking 
Af/y = Q- If we wish to account for all the metric fluctuations that are outside the horizon scale 
at a given time t, we should take Auv = o,{t)H, corresponding to a cutoff on physical momenta 
~ H. In this case, we find that the corrected formula for the path length takes the form 

A 



{S) = a{t)L 1 



367r 



a{t)HL + 



(6.16) 



For scalar fluctuations, A = iJ^/4e; similar contributions are obtained from tensor fluctuations, 
with A ~ H^. We find that the corrections to the path length due to the fluctuations rapidly 
become comparable to the length of the unperturbed path. For slow-roll inflation, this happens 
on a rapid time scale t ~ log[e/{H^L)], or when the satellite separation is of size 

a{t)L ~ eH-^ (6.17) 

(the result just accounting for tensor fluctuations is similar, without e). This effect completely 
overwhelms that of eq. (6.6), considered in [T7|[T8}|20]; the basic picture is that the curve can 
significantly shorten its length by taking advantage of the short-wavelength variations in the 
scale factor. While these effects might be removed by sufficient coarse-graining, taking the 
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scale q ~ we conclude that there are large fluctuations in the actual geodesic distances 
between satellites at these scales. These provide an apparent obstacle to a characterization of 
the geometry in terms of lengths of its geodesies, at distance scales longer than ~ And, 
if correlators are defined directly in terms of the proper distance in the three-geometry, as in 
(2.3) or eq. (6.8), this large fluctuating contribution likewise becomes significant P] 



6.2 Subleading effects 

This subsection argues that other contributions to the satellite separation are subleading to the 
effect in the preceding subsection; readers willing to take that for granted may wish to skip to 
the next section. There are multiple possible sources for such effects. The first category involves 
satellite dynamics: the possibility of an initial velocity for the satellites, or, in the quantum 
context, spreading of the wavepacket or even production of new satellites. The second category 
arises from metric fluctuations: satellites may move in response to metric perturbations, or 
there can be corrections to their separation from a one- loop contribution to {(). We consider 
these in turn. 



6.2.1 Satellite dynamics 

First consider the case where a satellite has an initial velocity in the comoving coordinates, 
V = dx^/dr. The Killing vector (2.4) then implies that 

k~ = a^{t)v (6.18) 
dr 

is conserved. Thus, any initial comoving velocity redshifts away as l/a^(t), and the satellite 
motion is rapidly dominated by the expansion. 

In a quantum framework, one should of course consider a wavepacket describing the satellite 
motion. Such a wavepacket is given by 

—f{k)u,{t)e''^- , (6.19) 

for some initial momentum-space profile /(k). Here the temporal wavefunctions are solutions 
of 

1 d I 3 dUk ( I / T 2 I 2 2\ 



a^— ^ ) + (r + mV)Mfc = . (6.20) 

0/ dib \ (Xib 

For dS, these take the form 

M,(t)=r73/2^p)(fcr/) (6.21) 



in terms of the conformal time (2.5), and with v"^ = 9/4 — m?/H'^. The positive frequency 
solution is given by H^'^\ On time scales short as compared to l/i?, we have \kri\ ^ 1, and an 



^^Onc might anticipate similar effects for tlie case of curves with constant acceleration, noted above, but 
checking this is left for future work. 
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initial wavepacket with width Ax ~ l/A/c <^1/H can be easily seen to evolve as in flat space, 
and in particular spreads as Ax{t) ~ tAk/m. At longer time scales, krj ^ 1, and one may use 
the asymptotics 



Hi'^ - -m ( f ) (6.22) 



2 Y 



to find 



71 \ki] J 

up to an overall constant. Thus, the spread Ax ~ 1/Ak in comoving coordinates x is frozen in 



time, This is consistent with the classical result that any initial velocity is rapidly redshifted 
to zero, having little subsequent effect. 

Finally, if satellites are treated as quantum objects, they may be created by fluctuations, 
just as with elementary particles, Boltzmann brains, or any other object in an inflating universe. 
If the number of microstates of such a satellite is TV, the production probability per unit four- 
volume is of size exp{— A/"}. Then, after N efolds, the total number produced is ~ exp{N —Af}. 
Since is bounded by the largest black hole in dS, and thus A^ <, S", we flnd a number of 
"Boltzmann" satellites comparable to the original number of satellites by ~ efolds. Before 
A^ ~ A/", this is a small effect. 



6.2.2 Force on satellites 



Satellites can also deviate from flxed comoving position due to forces exerted by the metric 
perturbations. To check that this effect is also small, flrst note that the satellite positions will 
be flxed in a different gauge, corresponding to geodesic-normal coordinates: 



AN = N' = . 



(6.24) 



Their motion in the comoving gauge (3.12) can then be found via the gauge transformation 

e lapse c 

A^^ = diijj 



relating these gauges. Specifically, in comoving gauge, the lapse and shift take the form [27] 

a 



N = 1 



with 




(6.25) 



(6.26) 



So, the gauge transformation (3.18) from (6.24) to comoving gauge is given by 



a : 



dt^C 



^^For H, the factor in front of the integral in (6.23) takes the form 7^^^ 
to Hubble dilution with massive oscillation. 



-3Ht/2-imt 



(6.27) 



corresponding 
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and 



C + — + h'^ / dt-djC (6.28) 



20,2 ad J CL 

and in the comoving gauge the satelhte trajectories are = —e*+ const. 



In (6.28), the last two terms redshift away. In slow roll, with 0o/^^ - /V"^ = 2e, the 



displacement is thus dominated by 



and for constant slow-roll parameter. 



= -e' ~ e^C , (6.29) 



d- 

x' = e-^C + const. (6.30) 



We can then evaluate the leading contribution to the displacement: 

-2 r „kfH^ r dk 



This is dominated in the infrared, leading to the dependence 

{{6x')^) ^ eH^L\ (6.32) 
and thus a typical small variation in separation given by 6L ~ y/tHL <^ L. 

6.2.3 The {C) tadpole contribution 

The tadpole of 7 must vanish, since it has to be invariant under 3-dimensional Euclidian 
transformations, and there is no such invariant, symmetric, traceless two-index tensor. However, 
one might still worry about possible IR divergent contributions from the tadpole of which in 



the slow-roll case could lead to effects competing with those in e.g. eq. (4.9). But if the splitting 



between background and fluctuations is properly handled, these effects will be subleading. This 



is most easily seen by first considering spatially- flat gauge, (3.13). In this gauge, the infiaton 



field can be written in terms of its backgound v.e.v. and the infiaton field fluctuations 

= 00 + (6.33) 

where by definition 

(0) = 00 , (^) = . (6.34) 



Now we can use the gauge transformation from spatially-fiat gauge into comoving gauge (3.12) 
on super-horizon scales to second order (see [27]) 



C = Cn + \{2e-r,)C (6.35) 
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with (n defined by the hnear relation 



H 

Cn = -^V^ • (6.36) 

00 



From the tadpole condition above, we then have (Cn) = 0, and to second order we obtain 

(C(x)) = ^(2e-r^)(C^(x)> . (6.37) 
Although the variance {("^{x)) contributes an IR divergence to the tadpole, it is clearly slow-roll 



suppressed compared to the other IR divergent contributions in eq. (4.9) 



6.2.4 Variance of the line integral 

Let us consider the contribution from tensor fluctuations to the variance of the line integral in 



eq. (4.8) around the torus in the slow-roll, since this is the most interesting case. Since the 
tadpole of the tensor modes vanishes, the variance is 

{ASAS) = ]a^it) lim // cix7...(t,x) / dx'^,,{t,x')) . (6.38) 
4 yz^y'z' \ Jq I 

When applying the boundary conditions for the torus, the momentum is discretized and we are 
projecting on the zero mode in the x-direction, obtaining 

(A5A5) = \a\t)L' ( -i I^^^l' ) ' (^.39) 

In de Sitter, we have |7fc(t)|^ ~ H'^/k^ for super-horizon modes, and the contribution to 
the variance becomes of order {AS AS) ~ a^L^H^, which is small and not IR enhanced. On 
the other hand, in slow- roll, there is a small IR enhancement, due to the non-fiat spectrum 
|7fc(t)|^ ~ Hl/k^ ~ {H'^ / k^){k / aH)^\ where nt = — 2e and e > 0. Here, we write the expansion 
rate as a function of momentum at horizon crossing instead of time by integrating eq.(5.18) 



of 21 



3 



8 In 



(6.40) 



where we assumed for simplicity a standard potential V{(j)) = A0^; "e" denotes values at the 



end of inflation. The sum ( 6.39[ ) is dominated by smallest k, and is of size 



(6-41) 

where c is an 0{1) constant. In order to evaluate the size of the effect, we assume that slow-roll 



inflation started in our local pocket near the self-reproduction scale, such that 21 



ffS 1 tt"^ , '7r2 

= IxN' = ^V'/' = — , (6.42) 
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where N is the total number of efolds, N^, ^ 60 is the number of e-folds left of inflation 
when the observable modes exits the horizon, and = (2.43 ± 0.11) x 10"^ is the observed 
amplitude of the primordial curvature perturbation spectrum. With these values one obtains 
(ifg/0g)iV^ ~ 10"^. Thus, the variance is of order a^L"^ x 10~^, which is still small. 



7 Volume fluctuations 

Another measure of spatial geometry in inflationary cosmologies is the volume; here we briefly 
comment on this as a diagnostic of fluctuations. 



First, in the flat gauge (3.13), the tracelessness of 'jij implies that the local volume element 



y/hd^x is just given by that in the unperturbed solution (3.3). Indeed, vanishing fluctuation in 



the local volume element can be thought of as a "physical" way of describing the origin of the 



choice of slicing (3.13). 



In contrast, when one chooses the time slicing specified by the comoving gauge (3.12) 



the fluctuations in the scalar field specifying the slicing lead to significant variations in the 



local spatial volume element. This effect has been explored in some detail in 22,23 , which, 
within the context of eternal inflation, specifically studied the probability for the volume on 
the reheating time slice, = (prh, to become infinite. Specifically, the fluctuation contribution 
to the volume element at quadratic order is given by 

(e3C)=et(^') , (7.1) 
and so the average volume of the reheating surface can be approximated by 

(V^(W)=e=^(^+^(^'(^)))(nt.n)) , (7.2) 
where trh is the time of reheating, tj„ is the initial time when inflation commenced, and is 



the number of efolds of the unperturbed solution (3.3). 

We can estimate when the effect on the volume at reheating becomes large; this happens 
when {("^{x)) becomes of order 2N/3. In a simple model of chaotic inflation with a monomial 
potential V{(j)) = the total number of e-folds scales with the initial inflaton field value, 
0i, as = 0?/(2a), while the variance scales as {("^{x)) = X(j)'^^^/{12Tr'^a^{a + 4)) |21j. This 



implies that the condition {(^{x)) > 2N/3 is satisfied when 

0,> (47rV(« + 4)/A)^ , (7.3) 
which is comparable to the inflaton field value at the end of the self-reproduction scale 

(psr ~ (127r2a7A)^ , (7.4) 
but smaller than the field value when the potential energy density becomes of order 0{Mp), 

0p ~ A-^/- . (7.5) 
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Thus, large fluctuations in the volume are closely associated with self-reproduction, between 
the Planck scale and self-reproduction scale, 0p > > (pgr, as described in 22,23 . In this 



regime, our perturbative approach breaks down 12 21 23 38 and there is no clear approach 
to computing the full "wavefunction of the universe." However, below this, in the non-eternal 



regime, the effect on the volume is small compared to the effect described in section 6.1 



8 Extensions and Conclusions 

In describing the quantum state of the geometry, one needs to formulate gauge-invariant q- 
observables; quantities such as {gfj,^{x)) or even (JZ{x)) are not gauge invariant. This paper has 
described one way of characterizing inflationary geometry, via nonlocal q-observables involving 
lengths of geodesies in the fluctuating spacetime. In the case where space is a torus, and in a 
time slicing with zero local deviation from classical expansion, we have found that fluctuations 
in the length around a cycle can become significant, and indeed such a curve can become 
shorter than an unperturbed curve, fixed in the original comoving coordinates, by a fractional 
amount that grows with the variance in the tensor fluctuations. This provides evidence that 
accumulation of these fluctuations is indeed an effect with physical consequences. Moreover, 
such contributions become large by a time scale t ~ (or t ~ RdsSds in terms of the de 

Sitter radius and entropy, in general dimensions). This is consistent with the suggestion that 
a perturbative description of dS fails at longer times, and with instability of dS space as a 



classical geometry 12 , 21 , 34 



We have thus found that one of the simplest formulations of gauge-invariant quantities 
receives large corrections due to quantum fluctuations. We might seek other characterizations 
of geometries, besides such geodesies, for better understanding of these and other effects. The 
challenge of respecting diffeomorphism invariance appears significant. For example, as noted, 
one may choose a path and examine the holonomy around that path, but the choice of path will 
not be diffeomorphism invariant. While workable measures of differences between geometries 
are scarce, we mention two possibly worth further exploration. The first is the Gromov-Hausdorf 
distance between two geometries. This is based on finding an embedding of these into a common 
geometry, and calculating a minimum distance over all such embeddings. A second is the 
DeWitt distance. If we have a fixed time slicing, so consider a difference 6g = gi ~ g2 between 
spatial metrics]^ this can be written 

\6g\^ = j (fx^g'^g^\5g,u5g,i + K5g,,5gu) ■ (8.1) 

Here is a constant; it may be chosen so that in a harmonic gauge, such as V^{6gij — ^gij6g) = 



0, the linear gauge transformation (3.6) has zero effect. DeWitt distance also appears to grow 
with the variance. 

If the inflationary spacetime is adorned with additional structure, say either quantum or 
classical matter, additional statements can be made. Specifically, matter allows one to formulate 



13 



One can also work directly in terms of the spacetime geometry. 
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other relational observables, as e.g. discussed in [2j. One may consider, as in the example of 
this paper, an initial distribution of (approximately) classical objects, and describe evolving 
features of the geometry with respect to these. 

In the case where we try to describe the geometry in terms of lengths between such "satel- 
lites," we find behavior that may have been unexpected. For starters, as the path endpoints 
reach separations large as compared to the Hubble scale 1/H, there is no spacetime geodesic, 
even in the unperturbed geometry, along which one can define the length. An alternative pre- 
scription is to define the length in terms of geodesies constrained to lie within spatial slices. 
Different choices of such spatial slices include those determined by either of the conditions 



(3.12) or (3.13). However, when one computes the lengths of the geodesies in the perturbed 
metric, one finds that the short-wavelength metric fluctuations lead to significant shortening 
of the geodesies, with an 0{1) contribution for geodesies of length S ~ 1/ (and, in the case 



of (3.12), shorter by a slow-roll factor e). Thus, one appears to lose control of a perturbative 
description of the geometry given in terms of such quantities]^ 

One may also try to formulate diffeomorphism-invariant correlators of local operators in 
terms of proper-distance separation, as in ( [23| |2||T6}]20)|24||25|. Here the above surprises also 
play a roll, forbidding such a construction with spacetime geodesies longer than ~ 1/if, and 
producing large fluctuations in geodesies constrained to spatial slices, when the operators are 
separated by distances 0{1/H^) (or RdsSds in general dimension). At longer length scales, 
characterization of correlators in terms of spatial geometry thus appears to encounter subtleties. 
While there are such nonlocal, diffeomorphism-invariant constructions that are sensitive to 
growth of fluctuations, there also do appear to be alternative sensible approximately local 
observables, that in particular are infrared safe, for observations restricted to small regions; 
speciflcally, those can be deflned with an effective resolution parameter that is the size of the 



horizon of a given observer, and apply to present-day observations [12]. (See also 13 14].) 

In short, when one attempts to describe the quantum state, and some of its observables, 
at sufficiently long scales, one flnds a breakdown in perturbative calculations. There is no 
known approach to accurately calculate this state, and this apparently requires nonperturbative 
gravitational dynamics]^ One can nonetheless apparently calculate certain physical quantities 
in small enough regions, over small enough time intervals. 
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^^This might be avoided at a coarse-grained level by introducing additional structure giving a long-distance 
averaging prescription, as noted above. 

A similar story appears present in black hole geometries, at times t ^ RS. This supports the resolution of 



15 



the information paradox proposed in 34 . This leaves an information problem, which requires understanding 



the long-distance modifications to the state resulting from the nonperturbative gravitational dynamics. 



29 



A Gauge fixing to higher order 

For our discussion of second-order effects in the H/Mp expansion, it is important to check gauge 
fixing to this order. Let us parameterize the metric as 

ds^ = -N^dt^ + hij{dx' + N'dt){dx^ + N^dt) ; hij = a\t){6ij + hij) . (A.l) 

Under combined time and space diffeomorphisms, 

t^t' = t-eO , ^^x" = x'-e' (A.2) 

the metric transforms to second order as 

Kj = hij + djei + ditj + 2t^H5ij + 82^^ , (A.3) 

where we expand N — 1, N'^ , hij, and e order-by-order in H/Mp, and 62hij is a term quadratic 
in the first-order fields and gauge parameters. 

For general hij (and lapse and shift), we want to find a gauge transformation so that 



1-3 



(A.4) 



with diiij = in = 0- 



From the trace of (A.3), we get to first order 

1 



(A.5) 



where H = d/a, and to second order 



,(2) 



1 

6H 



,(2) 



29,ef ^ + 62hi 



- (1)' (1)' 



(A.6) 



Now taking the divergence of eq. (A.3) gives to first order, using (A.5), 



(9% + ^9,9,)e« = -(a,/.S])-^9,/.;i)) 



(A.7) 



c.f (|3.9[)). This determines 'y^^^' to first order, via (|3.7[) and (|3.11[). 



At second order, we find 



1 



where J^j is given in terms of the second order h and the first order h and e's. This can 
be inverted, as in (3.11), to determine ej at second order. Note that this procedure may be 



iteratively continued to higher order. 
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